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We present measured J/tp production rates in d+Au collisions at ysjy„ — 200 GeV over a broad 
range of transverse momentum [pr ~ 0-14 GeV/c) and rapidity (—2.2 < y < 2.2). We construct 
the nuclear-modification factor RdAu for these kinematics and as a function of colUsion centrality 
(related to impact parameter for the -RdAu collision). We find that the modification is largest for 
collisions with small impact parameters, and observe a suppression (RdAu < 1) for pr < 4 GeV/c 
at positive rapidities. At negative rapidity we observe a suppression for pr < 2 GeV/c then an 
enhancement {RdAu > 1) for pr > "2 GeV/c. The observed enhancement at negative rapidity has 
implications for the observed modification in heavy- ion collisions at high pr- 

PACS numbers: 25.75. Dw 



I. INTRODUCTION 

Modifications of quarkonia yields when production 
takes place in a nuclear target, often termed cold-nuclear- 
matter (CNM) effects, give insight into the production 
and evolution of qq pairs. A number of effects are pre- 
dicted to occur in the presence of nuclear matter (for a 
recent review, see These include nuclear breakup, 

modification of the parton-distribution functions, initial- 
state parton-energy loss and, more recently, coherent 
gluon saturation. Measuring the production rate of 
quarkonia in a nuclear environment over a broad range 
of collision energies, and as a function of all kinematic 
variables, is the best way to disentangle these different 
mechanisms. 

The measurement of J/^ production rates over a broad 
range of rapidity (y) and transverse momentum (px) 
samples a wide range of parton momentum fraction (x) 
and energy transfer (Q^), providing a simultaneous con- 
straint on the modification of parton-distribution func- 
tions inside nuclei (nPDF's). The production of J/ip 
mesons, which at RHIC occurs mainly through gluon fu- 
sion, can provide critical input on the modification of the 
gluon distribution, which is probed only indirectly by the 
deep-inelastic scattering (DIS) data that forms the bulk 
of the current constraints on the nPDF paramctrizations. 

Measuring the px distribution of J/tp production al- 
lows access to pT-broadening effects, which are not con- 
strained by measurements of the rapidity dependence 
alone. The pT-broadening effects on quarkonia produc- 
tion at high energies are not well constrained by current 
data. New data for J/ip production over a broad range 
in pt is necessary to provide guidance for theoretical cal- 
culations. 

The CNM effects on J/tp production have been stud- 
ied in fixed-target p + A experiments at SPS, FNAL, and 
HERA @-[| spanning the center of mass energy range 
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y/s 17 — 42 GeV. The fixed-target results at midrapid- 
ity show greater suppression of J/ip production at lower 
collision energy [G;]. This has been interpreted as an 
increase of the nuclear breakup of the J/ ip through colli- 
sions with nuclei. At lower collision energy the crossing 
time of the nuclei is long enough for the J/ i/j to fully form. 
The fully formed J/ip has an increased probability of in- 
teracting with other nucleons in the collision, which can 
cause the breakup of the J/ip into heavy-meson pairs. At 
higher collision energies it is likely that the time required 
for the J/ip to fully evolve is as long, or longer than the 
crossing time of the collision. This may result in a de- 
crease in the probability of collisions with other nucleons, 
leading to less suppression of the J/ip production. 

The E866 and HERA-B Q experiments have mea- 
sured J/ip production as a function of px in fixed target 
p-\- A experiments. Results are presented in terms of the 
nuclear-suppression factor, a, which is obtained assum- 
ing that the cross section ioi p + A collisions scales as 
o'pA = o'pN X A"', where UpN is the proton- nucleon cross 
section and A is the mass number. They find a pt depen- 
dence of a, which is similar across a range of Feynman-x 
{xp) and Pt- At pt < 2 GeV/c they find a suppression 
in the J/ip production that transitions to an excess in 
the J/ip production at higher pT, which is characteris- 
tic of multiple scattering of the incident parton It 
is crucial to test these conclusions at the higher energies 
provided by d-l-Au collisions at RHIC in order to better 
understand the J/ip production mechanisms. 

Measuring, and understanding, the CNM effects on 
quarkonia production is critical to interpreting the results 
for J/ip production in nucleus-nucleus (A + A) colhsions. 
In 1986 Matsui and Satz predicted that the suppression 
of J/ip production in heavy-ion collisions would be a clear 
signature of the formation of a quark-gluon plasma [13] ■ 
The Debye color screening of the dense medium produced 
is expected to cause the dissociation of bound states, 
thereby causing a decrease in the observed production. 
Since then suppression of quarkonia production has been 
observed for a number of states, inclu ding th e J/ip and 
T, over a wide range in collision energy jTli - [l3 |. However, 
the interpretation of these results is still unclear. Before 
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the modification due to the produced medium can be de- 
termined, the CNM effects must first be corrected. This 
has been done at lower energies [ll|, but accurate data 
on CNM effects are still absent at the higher energies of 
RHIC and the LHC. 

Here we report new high-precision measurements of 
the J /ip production as a function of pr and collision cen- 
trality in d-|-Au collisions at ^s„„ = 200 GeV. We also 
present measurements of the J/V" -RdAu as a function of 
Pt, rapidity, and collision centrality using data for J /tp 
production in p+p collisions published in [l5j . PHENIX 
has previously measured the J/ip yield in d-|-Au colli- 
sions [IB, with data recorded in 2003. The data 
presented here, recorded in 2008, feature an increase in 
statistics of 30-50 times over those used in the previously 
published results, as well as a significant reduction of 
the systematic uncertainties. The rapidity dependence of 
J 1^ production in d-|-Au collisions from this data set has 
been previously published in 18]. This paper presents 
results for the px dependence of the J / t/j yield from the 
same data set. 



II. EXPERIMENTAL APPARATUS & DATA 
SETS 

The PHENIX detector [loj comprises three separate 
spectrometers in three pseudorapidity (77) ranges. Two 
central arms at midrapidity cover I77I < 0.35 and have an 
azimuthal coverage {(j)) of 7r/2 rad each, while muon arms 
at backward/forward rapidity cover —2.2 < 77 < —1.2 
(Au going direction) and 1.2 < ry < 2.4 {d going direc- 
tion), with full azimuthal coverage. 

In the central arms the J/ip yield is measured via 
diclcctron decays. Charged particle tracks are recon- 
structed using the drift chamber and pad chambers. Elec- 
tron candidates are selected by matching charged tracks 
to hits in the ring imaging Cerenkov (RICH) counters 
and clusters in the Electromagnetic Calorimeters (EM- 
Cal). In d-fAu collisions, a charged track is identified as 
an electron by requiring at least two matching RICH pho- 
totube hits within a radius of 3.4 < R[cm] < 8.4 with re- 
spect to the center defined by the track projection at the 
RICH. It is also required that the position of the EMCal 
cluster associated to the track projection match within 
±4(7, and that the ratio of the energy deposited in the 
EMCal cluster to the momentum of the tracks matches 
unity within ±2.5f7, where a characterizes the momen- 
tum dependent width of the matching distributions. A 
further cut of 200 MeV/c on the momentum of the elec- 
tron is added to reduce the combinatorial background, 
since the yield of electrons from J/tp decays observed in 
data and simulations is negligible below 200 MeV/c . 

At forward and backward rapidity, the J/ip yield is 
measured via dimuon decays. Muons are identified by 
matching tracks measured in cathode-strip chambers, re- 
ferred to as the muon tracker (MuTr) , to hits in alternat- 
ing planes of larocci tubes and steel absorbers, referred 



to as the muon identifier (MuID). Each muon arm is lo- 
cated behind a thick copper and iron absorber that is 
meant to stop most hadrons produced during the colli- 
sions, so that the detected muons must penetrate 8 to 11 
interaction lengths of material in total. 

Beam interactions are selected with a minimum-bias 
(MB) trigger requiring at least one hit in each of two 
beam-beam counters (BBCs) located at positive and neg- 
ative pseudorapidity 3 < |?7| < 3.9. The MB selection 
covers 88 ± 4% of the total d+Au inelastic cross section 
of 2260 mb [20]. 

The d-l-Au data sample used in this analysis requires 
the MB trigger to be in coincidence with an additional 
Level- 1 trigger. For electrons, this is a single electron 
EMCal RICH trigger (ERT), which requires a minimum 
energy deposited in any 2x2 group of EMCal towers, 
plus an associated hit in the RICH. Two thresholds on 
the minimum EMCal energy, 600 MeV and 800 MeV, 
were used, each for roughly half of the data sample. For 
muons, the level 1 trigger requires two tracks identified 
as muon candidates. The trigger logic for a muon candi- 
date requires a "road" of fired larocci tubes in at least 4 
planes, including the most downstream plane relative to 
the collision point. Additionally, collisions are required 
to be within ±30cm of the center of the interaction re- 
gion. Collisions in that range see the full geometric ac- 
ceptance of the central arms, and this cut also provides a 
reduction of the systematic uncertainties on the central- 
ity selection needed for the data from the muon arms. 
The data sets sampled via the Level- 1 triggers represent 
analyzed integrated luminosities of 62.7 nb~^ (electrons) 
and 54.0 nb~^ (muons) and nucleon-nucleon integrated 
luminosities of 24.7 pb^^ and 21.3 pb^^ respectively. 

TABLE L Characterization of the collision centrality for 
d+Au collisions along with the correction factor c (see text 



for details). 


Centrality 


(iVcoll) 


c 


c/(iVcoll) 


0-20 % 


15.1 ± 1.0 


0.94 ± 0.01 


0.062 ± 0.003 


20-40 % 


10.2 ± 0.7 


1.000 ± 0.006 


0.098 ± 0.004 


40-60 % 


6.6 ± 0.4 


1.03 ± 0.02 


0.157 ± 0.008 


60-88 % 


3.2 ± 0.2 


1.03 ± 0.06 


0.33 ± 0.02 


0-100 % 


7.6 ± 0.4 


0.889 ± 0.001 


0.117 ± 0.004 



The centrality, which is related to the impact parame- 
ter, b, of the d-t-Au collision is determined using the total 
charge deposited in the BBC located at negative rapidity 
(Au-going direction). The centrality is defined as a per- 
centage of the total charge distribution referenced to the 
greatest charge, i.e. 0-20% refers to the 20% of the total 
charge distribution with the greatest charge. On average 
the 0-20% centrality corresponds to collisions with the 
smallest b. 

For each centrality bin the mean number of nucleon- 
nucleon collisions ((iVcou)) is determined using a Glauber 
calculation f2H combined with a simulation of the BBC 
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d-Au at 200 GeV 

Centrality 00-20% 

Centrality 20-40% 

— Centrality 40-60% 
Centrality 60-88% 




III. J/iP ANALYSIS AND RESULTS IN THE 
MIDRAPIDITY REGION 



FIG. 1: (Color online) Nucleon-nucleon collision (A'coii) dis- 
tributions for each centrality bin obtained using a Glauber 
model for d+Au collisions described in the text. 



response (as described in [13]) • The resulting (A'^coii) val- 
ues for the centrality categories used in this analysis are 
shown in TablelH The A'coii distributions within each cen- 
trality bin are shown in Fig. [T] There is a significant 
overlap between the A'coii distributions for different cen- 
tralities. 

Also shown in Table H] is the correction factor c, which 
accounts for the correlation between the detection of a 
J / tp in the final state and an increase in the total charge 
collected in the BBC This correlation affects both 
the MB-trigger efficiency and the determination of the 
centrality of a given collision. The correction factors for 
each centrality bin are obtained within the same Glauber 
framework as the (A'coii) values by assuming that one of 
the N binary collisions produces a charge in the BBC 
that is characteristic of a hard-scattering process (the 
remaining N-1 binary collisions maintain a BBC charge 
distribution characteristic of soft scattering processes). 
The increase in the BBC charge from a hard process is 
tuned using real data. 

Since both c and (A'coii) are calculated in the same 
Glauber framework there are correlations between their 
uncertainties. These correlations are removed in the ratio 
of c/(A'coii), which occurs in the calculation of RdAn- The 
resulting values and uncertainties are given in the third 
column of Table HI The correction factor for 0-100% cen- 
trality contains an additional factor to extrapolate the 
measured yield, which covers only 88% of all d+Au col- 
lisions, to 100% of the d-l-Au inelastic cross section, es- 
sentially correcting for the efficiency of the BBC trigger. 
This correction is again determined within the Glauber 
framework using the parametrization of the BBC trigger 
efficiency. 



The procedure for analyzing the J — >■ e^e~ signal 
and the results in the central arms are discussed in this 
section. The extraction of the correlated e+e~ yield is 
discussed in Sec. IIII Al The estimation of the correlated 
background and losses due to the radiative tail in the J/ 
mass distribution is discussed in Sec. IIII Bl The estima- 
tion of the detector efficiencies is described in Sec. IIII CI 
The calculation of the J / ip invariant yield is detailed in 
Sec. IIIIDI The p+p baseline used in calculating RdAu is 
described in Sec. IIII El 



A. Correlated e^e Signal extraction 

The J/ip —> e~^e~ yield is measured using the invariant 
mass spectrum for all dielectron pairs where at least one 
of the electrons fired the ERT trigger. This selection is 
necessary to match the conditions under which the J/ip 
trigger efficiency is calculated (see Sec. IIII Cl) . An exam- 
ple of the dielectron mass spectrum is shown in Fig. [2] 
for 0-20% central collisions. In a given bin of px, ra- 
pidity, and collision centrality, the correlated e'^e~ yield 
(A'e+g- ) is determined by counting over a fixed mass win- 
dow of 2.8 < Mee [GeV/c^ ] < 3.3 the number of unlike- 
sign dielectrons, after the subtraction of the like-sign di- 
electrons, which arise by random association and so are 
representative of the combinatorial background within 
the unlike sign dielectron distribution. This method as- 
sumes that the acceptance is the same for e~ and e"*", 
which, while untrue at lower masses, is a good assump- 
tion in the J/tp mass range. At higher px, where statis- 
tical precision is limited, the yield, along with the sta- 
tistical uncertainties, are derived from Poisson statistics. 
Assuming both the unlike-sign (foreground) and like-sign 
(background) are independent, and assuming no negative 
signal, the combined distribution 



p(.) = 



A;=0 



fg 



fc)! 1 

bg\{fg-k)\ 2 



bg+fg-k ^k^-s 



k\ 



(1) 



represents the probability of a signal (s) given a number 
of unlike-sign dielectrons {fg) and a number of like-sign 
dielectrons {hg) (see [l5| for derivation). The mean and 
standard deviation of Eq. [T]are then used as the yield and 
uncertainties. A correlated e'^e~ yield in the mass win- 
dow 2.8 < Mee [GeV/c^ ] < 3.3 of approximately 8600 is 
obtained across all pt and collision centralities. 



B. Estimation of the Correlated Background and 
Losses Due to the Radiative Tail in the J/p Mass 
Distribution 

When using the like-sign subtraction method there re- 
mains a correlated background under the observed J /tp 
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FIG. 2: (Color online) (top) Invariant mass distribution of 
unlike-sign (filled circles) and like-sign (filled boxes) dielec- 
tron pairs in central d-f An collisions, integrated over pr and 
rapidity. Dashed vertical lines represent the mass range used 
to determine the correlated e^e~ yield, (bottom) Correlated 
dielectron invariant mass distribution for MB d-l-Au collisions. 
The line shapes are those used to extract the continuum and 
radiative tail contributions to the correlated e'^e~ yield in the 
mass range 2.8 < Me, [GeV/c^ ] < 3.3 . 



peak. This background comes mainly from open-heavy- 
flavor decays and Drell-Yan pairs, and must be separated 
from the J/ifj signal of interest. Counting the dielectron 
signal only over a fixed mass window also causes an un- 
derestimate of the J/ip yield due to the fraction of the 
J/ip line shape that falls outside the mass window of 
choice. These two effects are quantified by using simu- 
lated particle line shapes fitted to the real data distribu- 
tion. 

The J/tp and ^' mesons with uniform distributions in 
PT (0 < pT [GeV/c ] < 12) and rapidity {\y\ < 0.35) are 
decayed to e~^e~ and the external radiation effects are 
evaluated using a GEANT-3 based model of the PHENIX 



detector (described in Sec. lIII Cl) . While a uniform distri- 
bution in pt is unrealistic, the J/ip rapidity distribution 
is roughly constant within \y\ < 0.35. When used here, 
the J/ip and ^' line shapes will be compared to px inte- 
grated data as a function of invariant mass only, with a 
mass resolution fitted to the data, and therefore the ef- 
fect of using a uniform px distribution is negligible. The 
line shape for J/4> radiative decays {J/tp e'''e~7), also 
called internal radiation, is based on calculations of the 
mass distribution from QED [22*1 convoluted with the de- 
tector resolution. 

Line shapes for the correlated background from heavy- 
flavor decays along with Drell-Yan pairs are simulated us- 
ing PYTHIA [2^ . The correlated background from heavy- 
flavor decays comes from semi-leptonic decays of corre- 
lated DD and BB (i.e. D ^ e+ + X and D ^ e- + X). 
The decay electrons from pythia are then run through 
the same GEANT simulation of the PHENIX detector to 
evaluate the external radiation effects. These line shapes 
are generated assuming p + p collisions, and no correc- 
tions for CNM effects (i.e. application of nPDF modifica- 
tions) are applied to the distributions. We assume that 
the CNM effects on these distributions are likely small 
and roughly constant over the narrow mass window used 
due to the x values probed. No suppression of heavy- 
flavor production has been observed in d-|-Au collisions, 
and we assume that any suppression, if it exists, does not 
significantly effect the overall line shapes. 

The line shapes are then fitted to the px and colli- 
sion centrality integrated invariant mass spectrum over 
the mass range 2 < M^e [GeV/c^ ] < 8 where the nor- 
malizations on the J/4', "0'; correlated heavy flavor, and 
DY are free to vary independently. The best flt is shown 
in the quarkonium mass region in Fig. [21 where the con- 
tinuum line shape is the combination of correlated e~'"e~ 
pairs from DD, BB, and DY decays, and the J/tp and tp' 
line shapes are the combinations of the line shapes from 
both internal and external radiation effects. Within the 
mass window 2.8 < Al^e [GeV/c^ ] < 3.3 the correlated 
continuum contribution (econt) is found to be 6.6 ± 0.2% 
and the fraction of the J/^ line shape contained within 
the mass window (crad) to be 94.3 ± 0.2%, where the 
uncertainties are derived from the uncertainty in the fit. 
The disagreement between the fit and the data in the 
3.7 < Mee [GeV/c^ ] < 4.5 mass range is likely due to 
the inability of the DD and BB line shapes to match the 
shape of the data at higher mass. However, large changes 
in the ratio of their contributions have only a small ef- 
fect on the extracted values of Ccont and Siad , and this is 
accounted for in the quoted uncertainties. 



C. Acceptance and Efficiency Studies 

The J/ip acceptance is investigated using a GEANT- 
3 [13 based Monte Carlo model of the PHENIX detector. 
Dead and malfunctioning channels in the detector are re- 
moved from both the detector simulation and real data. 
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The accuracy of the simulations is tested by comparing 
simulated single electron distributions with those from 
real data. The agreement across the detector and data 
taking period is determined to be within 3.2%. A con- 
servative estimate, which assumes that the uncertainty is 
correlated for both electrons in a pair, of 2 x 3.2% — 6.4% 
is assigned as a systematic uncertainty on the J/ip ac- 
ceptance based on the quality of the matching between 
simulations and data. 

To determine the J/ip acceptance, J/ip — )■ e+e" decays 
are simulated with uniform distributions in px, rapidity 
(|y| < 0.5) and collision vertex. While distributions uni- 
form in pt are not realistic, the corrections are made 
over a small pT bin where the real distribution can be 
approximated as linear. This assumption, and the effect 
of bin sharing, is tested later and taken into account when 
assigning systematic uncertainties. The fraction of J ftp 
decays that are reconstructed corresponds to the combi- 
nation of the geometric acceptance and the electron ID 
efficiency {A x EcId)- The resulting A x EoID is shown as a 
function oi pr in Fig. [31 It has an average value of 1.5% 
in 1 unit of rapidity. The dip in A x eeiD followed by a 
continual increase with pT marks the transition from the 
e~^e~ pair at low pT being produced back to back and be- 
ing detected one in each of the PHENIX central arms, to 
the pair at high pT being produced in a collinear manner 
and being detected both in the same PHENIX central 
arm. The low point at pt ~ 3 GeV/c corresponds to 
the e+e~ being produced at roughly 90° relative to each 
other in the lab frame, which due to the PHENIX geome- 
try has the lowest probability for detection. The electron 
ID efficiency, which is mainly due to track reconstruc- 
tion cuts used to avoid the misidentification of hadrons 
as electrons, was cross checked using electrons from 7r° 
Dalitz decays and 7 conversions as described in [isj , and 
a systematic uncertainty of 1.1% is assigned based on 
that comparison. The effect of momentum smearing on 
the electrons in simulations, which can cause a J/ip to be 
reconstructed into a different pt bin than the one it was 
generated in, was also investigated. The effect was found 
to be minimal for all but the highest px bins and an un- 
certainty of 0.2% was assigned based on a Monte-Carlo 
study effect and a parametrization of the measured mo- 
mentum resolution for electrons. A combined uncertainty 
of 6.5% is assigned to the J/4> A x eeio by adding the 
simulation/data matching, elD, and momentum smear- 
ing uncertainties in quadrature. 

The ERT trigger efficiency is evaluated using simula- 
tions of J/ip decays and parametrizations of the single 
electron trigger efficiencies in each trigger tile. A MB 
data sample of single electrons is used to measure the pt 
dependent efficiency of each 2x2 EMCal trigger tile and 
each RICH trigger tile independently by calculating the 
fraction of electrons that fired the trigger tile compared 
to all those passing through it. The resulting distribu- 
tions are then fitted with an error (uniform) function for 
each trigger tile in the EMCal(RICH). These functions 
are then used with simulated J/iJj decays to estimate the 



efficiency of the ERT trigger for triggering on e~^e pairs 
from J ftp decays (e^^^). The trigger efficiency is evalu- 
ated only for simulated J/ip decays for which both elec- 
trons passed an acceptance and trigger check in order 
to avoid double counting the acceptance efficiency. This 
procedure is repeated independently for each of the two 
EMCal trigger thresholds used during the run. The pr 
dependence of e'^^ is shown in Fig. [3l where both ERT 
trigger thresholds have been combined using the relative 
luminosities of each data sample. It has an average value 
of 77%. The dip seen at pr ~ 3 GeV/c is due to the kine- 
matics of the J/ip decays. In that px range there is a high 
probability for the decay electrons to have unbalanced 
momenta, where one of the electrons will have a momen- 
tum below or near the trigger threshold, resulting in a 
lower probability for triggering on the J/ip- The effect 
of the fit function used in the EMCal trigger tile efficien- 
cies is investigated by replacing the error function with a 
double-Fermi function. This gives an average change in 
the J/ip ERT efficiency of 0.31%. The statistical uncer- 
tainty in the trigger tile efficiency leads to an uncertainty 
in the J/ip ERT efficiency of 1.6%. Summing these uncer- 
tainties in quadrature gives a total uncertainty on 
of 1.6%, which is heavily dominated by the uncertainty 
in the efficiency of each ERT trigger tile. 

The detector occupancy effect is negligible, even in 0- 
20% central d-l-Au collisions (a finding consistent with 
previous embedding studies in peripheral Cu-|-Cu [2^ 
with similar multiplicities). A 1% systematic uncertainty 
was assigned based on studies where simulated J/ip de- 
cays were embedded into real events. This result agrees 
well with the studies done in [l^ , where a slightly larger 
systematic uncertainty was assigned because of the lower 
statistical precision of the simulations used. 
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FIG. 3: The J/ip acceptance x electron ID efficiency (a) 
and J/^p ERT trigger efficiency (b) as a function of pr for 
\y\ < 0.35 where the shaded boxes represent the systematic 
uncertainties. 
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D. Invariant Yield Results 

The J 1^ invariant yield in a given rapidity, transverse- 
momentum, and centrality bin is 



1 



cN 



Bii d?N 
2ttpt dydpT 2'kpt^Pt^V ^EVTCtot ' 



(2) 



where Ba is the J/-!/; — ?> l^l^ branching ratio, Nj/.^ is 
the measured J /ip yield, A'evt is the number of sam- 
pled MB events in the given centrality bin. Ay is the 
width of the rapidity bin, Ap^ is the width of the pt 

bin, etot = ^ x CeiD eEprr ^rad and c is the BBC bias 
correction factor described in Sec. At midrapidity 
Nj/^ = iVg+e- (1 — econt)j where Ccont is the correlated di- 
electron continuum contribution in the J mass range. 
The 0-100% centrality integrated J/tp invariant yield is 
shown as a function of in Fig.|4l and for four centrality 
bins in Fig. [5j Here the values shown represent the aver- 
age over the pr bin and are plotted at the center of the 
bin, as this provides the measured information without 
introducing further systematic uncertainties. 



TABLE IL A summary of the systematic uncertainties at 
\y\ < 0.35. 



Source 



Value 



Type 



Embedding 

trad 
Econt 

c (Invariant yield only) 
c/ (Afcoii) (-RdAu only) 

A X EcID. 

J/i' 

•^ERT- 

Stat. Uncertainty on 
the correlated e^e~ yield 



1.0% 
0.2% 
0.2% 
0.1-5.8? 
3-6% 
6.5% 
1.6% 
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FIG. 4: (Color Online) J/il> invariant yield as a function 
of Pt for p+p and 0-100% centrality integrated d+Au col- 
lisions. The type C systematic uncertainty for each distribu- 
tion is given as a percentage in the legend. The midrapidity 
d+Au and p+p results are discussed in Sees. IIIIDI & IIIIEI 
while the forward/backward rapidity results are discussed in 
Secs. UVCl fc UVDl 



to be zero. To remain consistent, this effect is removed 
from the p+p baseline as well, so that, assuming the po- 
larization does not change drastically between p+p and 
d-|-Au, the effects will cancel in the nuclear ratio, RdAu- 
The p+p invariant yields as a function of pr used in this 
work, shown in Fig. 21 have been converted from the in- 
variant cross sections published in [l^ using an inelastic 
cross section of 42 mb. 



A summary of all the relevant systematic uncertain- 
ties at midrapidity is shown in Table |TT1 along with their 
classification into Type A, B, or C uncertainties. Type A 
represents uncertainties that are uncorrelated from point 
to point. Type B represents uncertainties that are corre- 
lated from point to point, and Type C represents uncer- 
tainties in the overall normalization. 



E. p+p Baseline 

The p+p baseline used to calculate -R^au is extracted 
from 2006 data published in [l^. The integrated lumi- 
nosity was 6.2±0.6 pb^^. In the analysis, described in 
detail in [l^, the effect of the J/ip polarization on the 
J/ip acceptance is included. This effect is not included in 
the d-|-Au result presented here due to a lack of knowl- 
edge of the effects of a nuclear target on the J/tp po- 
larization. The J/tp polarization is therefore assumed 



IV. J/iP ANALYSIS AND RESULTS IN THE 
FORWARD/BACKWARD RAPIDITY REGION 

The procedure for analyzing the J/^p — >■ A*"*"/^" sig- 
nal at backward and forward rapidity in the muon arms 
is discussed in this section. The procedures are simi- 
lar to those detailed in jTsj . with only a brief summary 
presented here, except where there are differences. As 
in [isj , the rapidity region of the forward muon arm used 
in the analysis was truncated to 1.2 < y < 2.2 to match 
the rapidity coverage of the backward muon arm. The 
extraction of the raw J ftp yield is discussed in Sec. lIVAl 
The estimation of the detector efficiencies is described in 
Sec. IIVBI The calculation of the J/ip invariant yield is 
detailed in Sec. IIV CI The p+p baseline used in calculat- 
ing RdAu is described in Sec lIVDl 
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A function consisting of an exponential component 
combined with two Gaussian distributions, which are 
used to better reproduce the mass resolution present in 
the muon arms, was used to fit the dimuon mass dis- 
tribution, convolved with a function to account for the 
variation in acceptance over the invariant mass range. 
An example of the fitted mass distribution is shown in 
Fig. ini Both the J/V' component of the fit, and direct 
counting after the subtraction of the fitted exponential 
background, are used to evaluate the J/^ yield. The 
difference between the two methods is taken as a Type 
A systematic uncertainty. This uncertainty is typically 
small (« 2%) but can be significantly larger at high pT 
where there are fewer counts. Measured J/tp yields of ap- 
proximately 38000 and 42000 are obtained at backward 
and forward rapidity, respectively. 



FIG. 5: (Color Online) J/^ invariant yield as a function of 
Pt for each centrality at \y\ < 0.35. The type C systematic 
uncertainty for each distribution is given as a percentage in 
the legend. 



A. J/ip — > /i^/i Signal extraction 

At forward and backward rapidity, the invariant mass 
distribution is calculated for all unlike-sign dimuons in 
events that pass the trigger requirements described in 
Secini The combinatorial background is estimated from 
the invariant mass distribution formed by pairing unlike- 
sign muon candidates from different events. This is done 
to reduce the background statistical uncertainty below 
what is possible by subtracting like sign pairs from the 
same event, and is needed because the signal to back- 
ground present at forward/backward rapidity is smaller 
than at midrapidity. The mixed event muon pairs are 
required to have vertices that differ by no more than 3 
cm in the beam direction. The mixed event spectrum is 
normalized by the factor 



rmxed\ 
11- 11- 



/ f ATvaixed\ I ATmixed\ 



(3) 



where 7V"™'= and iV""'= 



are the number of pairs formed 
from two muons in the same or in mixed events, respec- 
tively. 

The remaining correlated dimuon mass distribution 
after the subtraction of the mixed event combinatorial 
background contains dimuons from J/ij} and '0' decays, 
as well as correlated dimuons from heavy-fiavor decays 
and Drell-Yan pairs. Due to the momentum resolution 
of the detector, there is no clean discrimination between 
the J/ij} and V'' in the mass distribution. However the 
■0' contribution is expected to be negligible in the mass 
window of interest. 
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FIG. 6: (Color Online) (top) Invariant mass distribution of 
unlike-sign (filled circles) and like-sign (filled boxes) dimuon 
pairs for 2 < pr < 2.25 GeV/c at forward rapidity and 60- 
88% central events, (bottom) Invariant mass distribution of 
correlated dimuon pairs after the subtraction of the com- 
binatorial background. The solid line represents the fit to 
the invariant mass distribution, which includes the double 
Gaussian signal component (dot-dashed line) and exponen- 
tial background (dotted line). 
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B. Acceptance and Efficiency Studies 

Studies of the response of the muon arm spectrome- 
ters to dimuons from J/'0 decays are performed using a 
tuned GEANT3-based simulation of the muon arms, cou- 
pled with a MuID trigger emulator. The MuID panel- 
by-panel efficiencies are estimated using the fraction of 
reconstructed roads in real data. Where statistics are 
limited, the operational history of each channel recorded 
during the run was used to estimate the efficiency. A 
systematic uncertainty of 4% is assigned to the MUID ef- 
ficiency based on this comparison. Charge distributions 
in each part of the MuTr observed in real data, along 
with dead channels and their variation with time over 
the run, are used to give an accurate description of the 
MuTr efficiency within the detector simulation. 

The J/'0 acceptance x efficiency {Axe) evaluation uses 
a PYTHIA simulation with several parton distributions as 
input to account for the unknown underlying rapidity de- 
pendence of the J/ 4> yield. A 4% systematic uncertainty 
is assigned based on changes in the input parton distri- 
butions. A systematic uncertainty of 6.4(7)% on the J/ip 
yield is assigned to the backward(forward) rapidity due 
to the uncertainties in the acceptance x efficiency deter- 
mination method itself. 



TABLE III: 
1.2 < \y\ < 2.2. 



Invariant Yield Results 



The dominant systematic errors at 



Source 


Value (S/N) 


Type 


c (Invariant yield only) 


0.1-5.8% 


C 


c/ (Acoii) {RdAvL only) 


3-6% 


C 


MC Input Distributions. 


4% 


B 


MuTr Efficiency 


2% 


B 


MUID Efficiency 


4% 


B 


Acceptance 


6.4/7% 


B 


Fit Type 


«2% 


A 


Stat. Uncertainty on 
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FIG. 7: (Color Online) J/ip invariant yield as a function 
of pt for each centrality for a) —2.2 < y < —1.2 and b) 
1.2 < y < 2.2. The type C systematic uncertainty for each 
distribution is given as a percentage in the legend. 



D. p+p Baseline 



The J/iJj invariant yield at backward/forward rapidity 
is calculated using Eq. [5J where etot — Axe. A summary 
of the systematic uncertainties is given in Table IIIIl The 
backward and forward 0-100% centrality-integrated J/ip 
invariant yields are shown as a function of pt in Fig. 21 
while the J/^ invariant yields are shown as a function of 
Pt in each centrality bin in Fig. [71 



The p+p baseline used to calculate i?dAu is extracted 
from a combined analysis of data taken in 2006 and 2008, 
published in [Tsj . The combined integrated luminosity 
was 9.3±0.9 pb~^. As discussed in Sec. IIIIE"! the effect 
of the J/ ijj polarization on the J/ ip acceptance is removed 
from the results used here. The J/ip invariant yield in 
p+p collisions at forward/backward rapidity used here 
is shown as a function of px in Fig. 21 for convenience, 
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where we have converted from the invariant cross sections 
pubhshed in [l^ using an inelastic cross section of 42 mb. 



V. CALCULATION OF (P|) 

The Ijr^) is calculated for each of the J/i/' invariant 
yields presented in Sees. IIII Dl and lTVCl and the resulting 
values are shown in Table HVl 



function of the form 



TABLE IV: (py) results for p+p and d+Au collisions where 
the first quoted uncertainty corresponds to the type A uncer- 
tainties and the second corresponds to the type B uncertain- 
ties. 



System 



y range 



Centrality Ij)^) [GeVVc^] 



P+P 
P+P 

d+Au 
d-l-Au 
d+Au 

d4-Au 
d-l-Au 
d-l-Au 
d-l-Au 

d-l-Au 
d-l-Au 
d+Au 
d+Au 

d+Au 
d+Au 
d+Au 
d+Au 



1.2 < \y\ < 2.2 
\y\ < 0.35 

-2.2 <y < -1.2 
\y\ < 0.35 
1.2 <y < 2.2 

-2.2 < y < -1.2 
-2.2 <y < -1.2 
-2.2 <y < -1.2 
-2.2 <y < -1.2 

\y\ < 0.35 
\y\ < 0.35 
\y\ < 0.35 
\y\ < 0.35 

1.2 < y < 2.2 

1.2 < y < 2.2 

1.2 < y < 2.2 

1.2 < y < 2.2 



0-100% 
0-100% 
0-100% 

0-20% 
20-40% 
40-60% 
60-88% 

0-20% 
20-40% 
40-60% 
60-88% 

0-20% 
20-40% 
40-60% 
60-88% 



3.64+0.03+0.06 
4.46+0.14+0.18 

4.09+0.06+0.09 
5.101E;;1o±0.11 
4.05+0.05+0.10 

4.22+0.08+0.09 
4.06+0.08+0.09 
4.01+0.09+0.09 
3.92+0.10+0.09 



5.24 
5.27 



-1-0.19 
-0.16 
-1-0.22 



+0.10 
+0.12 



3-)-0.29_ 
'-0.26- 

4.60lnl4±0.15 



4.15+0.06+0.10 
4.13+0.07+0.11 
3.94+0.07+0.10 
3.80+0.08+0.10 



Unlike in previous analyses [16|, where the (p^) was 
calculated for < 5 GeV/cdue to statistical limita- 
tions at high pt, here we have calculated the {p'^) over 
the full Pt range. First the (p^) was calculated numer- 
ically up to the Pt limits of the measured distribution 
((Pt) lpT<p"'"') • The correlated uncertainty was prop- 
agated to (Pt) |pT<p"'"' by sampling the type B uncer- 
tainty distributions of the first and last pt point of the 
invariant yield, and assuming a linear correlation in be- 
tween. For a more detailed description of this procedure 
see Appendix I A 31 

To account for the differences in the pt limits of the 
various distributions, the {Pt)\pt<p^'''^ value was cor- 
rected to the pt range from zero to infinity. This was 
done by fitting the distribution with a modified Kaplan 



Pt 
Pi. 



2\ P2 



where each parameter was free to vary. The ratio 



(p|>[0,pt'""''] 



(4) 



(5) 



was then calculated from the fit and applied to the nu- 
merically calculated {pt)\pt<p^'^'^- all cases the cor- 
rection factor was small (fc < 1.03), and an uncertainty 
in the correction factor based on the fit uncertainty is 
included in the Type B uncertainties shown in Table IIVI 
For a more detailed description of this procedure, includ- 
ing the fit results and the calculated values of k see Ap- 
pendix |X1 

The (pt) for p+p collisions was previously published in 
[isj . But we report the result here with the effect of the 
J/tp polarization on the acceptance removed. The results 
are in good agreement with those presented in and 
are shown in Table IIVI 

Figure [8] shows A(p|,) = {pT)dAu ~ iPT)pp as a function 
of A^coU- There is a broadening in the pt distribution 
with respect to p+p, which increases with iVcou, and is 
similar at forward and backward rapidities. We observe 
a larger increase in the pt broadening at midrapidity. 
However, this observation is tempered by the relatively 
large uncertainties present in the data. 

1-2| 1 I I I I I I I I I 1 1 I I I I I I I 1 1 I 1 I I I I I I 
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FIG. 8: (Color Online) The difference between the .7/tp (pt) 
in d+Au and p+p collisions as a function of A^coii in d+Au 
collisions. The boxes drawn at A{py) = represent the com- 
bined statistical and systematic uncertainties from the p+p 
calculation. 



VL THE J/t/) RdAu 

To quantify the d+Au cold nuclear matter effects, the 
J/f/' RdAu is calculated for a given pt, y, and centrality 
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FIG. 9: (Color Online) J/tp nuclear modification factor, RdAu, 
as a function of pr for (a) backward rapidity, (b) midrapidity, 
and (c) forward rapidity 0-100% centrality integrated d+Au 
collisions. Curves are model calculations detailed in Sec. I VI Al 



bin as: 



e d^N:}l^^\t)/dydpT 
(A^con(*)> dmP+P/dydpT ' 



(6) 



where d^ N'^'^^^{i) / dydpr is the d+Au invariant yield for 
the i*^ centrality bin, d'^N^'^^ / dydpr is the p+p invariant 
yield for the same pr and y bin, and {Nco\i{i)) is the av- 
erage number of binary collisions for the given centrality 
bin, as listed in Table ID 

The 0-100% centrality integrated J/ip RdAu as a func- 
tion of pt is shown in Fig. [S] for each of the three ra- 
pidity regions. The numerical values can be found in 
Table I VIII IVIIIl and IIXI for backward, mid and forward 
rapidity, respectively. Figure [9] shows a different behav- 
ior for RdAu at backward (—2.2 < y < —1.2) as opposed 
to mid {\y\ < 0.35) and forward (1.2 < y < 2.2) rapidi- 
ties. At backward rapidity, the -RdAu is suppressed only 
at the lowest px, with a rapid increase to RdAu — 1-0 at 
Pt ~ 1.5 GeV/c . The mid and forward rapidity data, 
on the other hand, exhibit a similar level of suppression 
at the lowest px, but a much more gradual increase in 



-RdAu with Pt, increasing to i?dAu = 1.0 only at ~ 4.0 
GeV/c . Figure [TU] shows the same 0-100% RdAu vs pr 
at all rapidities overlaid. It is striking that the shape and 
absolute scale for the mid and forward rapidity data is 
nearly consistent across the entire pt range of the data. 
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FIG. 10: (Color Online) J/ip RdAu, as a function of pt for 0- 
100% centrality integrated d-|-Au collisions at each rapidity. 
The Type C systematic uncertainty for each distribution is 
given as a percentage in the legend. 

Due to the statistical limitations of the data at high 
Pt, it is unclear from Fig. [9] whether the i?dAu increases 
significantly above one. To investigate the high-pT be- 
havior of the i?dAu at each rapidity, the average RdAu 
was calculated for > 4 GeV/c by fitting each distri- 
bution with a constant. The results are shown in Table IVl 
along with the fit uncertainties, which take into account 
only the type A uncertainties on the data. Since the type 
B uncertainties are roughly consistent in the fit range, we 
have chosen here to add the average type B uncertainty 
for Pt > 4, GeV/c in quadrature with the type C un- 
certainty. We find that at mid and forward rapidity the 
average RdAu for pt > 4 GeV/c is consistent with 1.0, 
while at backward rapidity the average i?dAu is greater 
than 1.0. 

The production of a J/ip at forward rapidity in ^ -I- A 
collisions involves a low- a; gluon colliding with a high- 
X gluon. The symmetry due to identical colliding nuclei 
results, essentially, in the folding of the forward and back- 
ward rapidity RdAu- The production of a J/tp at midra- 
pidity results, essentially, in the folding of the midrapid- 
ity RdAu with itself. This picture is simplistic and leaves 
out many details, but it gives some expectation for the 
result of the modification of J/ip production in A + A 
collisions due to CNM effects. If we extrapolate the ob- 
served behavior of RdAu to the modification of J ftp's pro- 
duced at forward rapidity in A + A collisions, we would 
expect a i?AA contribution from CNM effects to be simi- 
lar to, or greater than, 1.0 at high pT and a modification 
similar to 1.0 at midrapidity. The observation at midra- 
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pidity of a J /ip Raa in Cu+Cu collisions that exceeds, 
but is consistent with, 1.0 at high px [1^ may there- 
fore be largely accounted for by the contribution from 
CNM effects. Further work is needed to understand the 
detailed propagation of measured results in d+Au colli- 
sions to an expected CNM contribution in A+A collisions 
before this can be fully understood. 



TABLE V: The average 0-100% RdAu for pr > 4 
GeV/c where the first quoted uncertainty corresponds to the 
fit uncertainty and the second corresponds to the combined 
type B and C systematic uncertainties. 



Rapidity 


i?dAu (PT > 4 GcV/c ) 


-2.2 <y < -1.2 


1.27±0.06±0.11 


\y\ < 0.35 


0.97±0.14±0.16 


1.2 < y < 2.2 


1.03±0.06±0.11 



Figures [TTirT^ and [T^ show Rd Au vs pr in four central- 
ity bins for backward rapidity, midrapidity, and forward 
rapidity, respectively. Numerical values can be found in 
Tables Kl El Kn] and Km] for 0-20%, 20-40%, 40-60% 
and 60-88% central collisions respectively. For peripheral 
collisions the RdAu remains consistent with 1.0 within 
statistical and systematic uncertainties across all pr in 
all rapidity regions. 



A. Comparison with Model Predictions 

As mentioned previously, various models have been 
suggested to describe the cold nuclear matter effects on 
J /ip production. The models that will be discussed here 
include a combination of physical effects such as shadow- 
ing, nuclear breakup, and the Cronin effect. 

Shadowing, the modification of the parton distribu- 
tions within a nucleus, is calculated using parametriza- 
tions of deep inelastic scattering data in the form of nu- 
clear modified parton distribution functions (nPDF's). 
There are a number of nPDF sets available, includ- 
ing nDSg 28], EKS98 [29] and EPS09 [30], which pro- 
vide distributions of this modification based on differ- 
ent parametrizations of the available data. For J /^p pro- 
duction in d-l-Au collisions the relevant distributions are 
those providing the modification of the gluon distribu- 
tion within a Au nucleus, as J/ip^s are produced pri- 
marily through gluon fusion at y^s^ =200 GeV. The 
nPDF's provide modifications as a function of parton mo- 
mentum fraction [x) and energy transfer (Q^). Knowl- 
edge of the J/V' production kinematics is then needed 
to produce a modification to J /ip production in d-|-Au 
collisions. For J /ip production at backward rapidity and 
< pt < 8 GeV/c , a range of roughly 0.051 < a; < 0.39 
in the Au nucleus is probed, assuming simple 2—^1 kine- 
matics. While 2— >-l kinematics are inadequate to de- 
scribe the production of a J/ip with nonzero pr, they 
are used here to provide a simple estimation of the x 
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FIG. 11: (Color Onfine) J /tp — >■ /i^/i RdAu, as a function of 
Pt for a) central, b) midcentral, c) midperipheral, and d) pe- 
ripheral events at -2.2 <y < -1.2. The 60-88% RdAu point 
at Pt = 5.75 GeV/c has been left off the plot, as it is above 
the plotted range and has very large uncertainties, however it 
is included in Table IXlIIl Curves are calculations by Lansberg 
et al. [231 discussed in the text. 



and ranges covered. Likewise, midrapidity covers 
roughly 0.0094 < a; < 0.071 and forward rapidity cov- 
ers roughly 0.0017 < a; < 0.013. A range of roughly 
10 < (3^[GeV^/c^] < 74 is probed at each rapidity under 
the same assumptions. The data thus provide a strong 
constraint to shadowing models over a wide range of x 
and Q^. 

Nuclear breakup, the dissociation of cc pairs that 
would have formed J/tp^s through collisions with nucle- 
ons, is often parametrized through a breakup cross sec- 
tion. Little theoretical or experimental guidance cur- 
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FIG. 12: (Color Online) J/ip e^e~ RdAu, as a function 
of pt for a) central, b) midcentral, c) midperipheral, and d) 
peripheral events at \y\ < 0.35. Curves are calculations by 
Lansberg et al. [2^ discussed in the text. 



FIG. 13: (Color Online) J/tp ^ fJ-^ 1^^ RdAu, as a function 
of Pt for a) central, b) midcentral, c) midperipheral, and d) 
peripheral events at 1.2 < y < 2.2. Curves are calculations by 
Lansberg et al. [2^ discussed in the text. 



rently exists on the exact nature of this effect due to 
the many complications and competing effects involved 
in J/ ip production in p{d)+A collisions. Often this effect 
is modeled by a simple "effective" cross section, vi^hich re- 
mains constant with px, however there are a number of 
models, including a dynamic breakup cross section that 
changes based on the kinematics of the produced J/ip. 

The broadening of the px distribution, termed the 
Cronin effect [31], is typically attributed to multiple elas- 
tic scattering of the incoming parton before the hard col- 
lision that produces the J/tp- This modifies the pr de- 
pendence of the J/ip production by adding pT vectorially 
to the incoming parton. This generally causes a decrease 
in J/ijj production at low px and a compensating increase 



at higher pT [pt ~ 5 — 10 GeV/c ), which eventually falls 
off at yet higher px [pt ~ 10 GeV/c ). 

The first set of model calculations that we discuss is 
by Kopeliovich et al. jH, calculates the effects on 
a cc dipole propagating through a nucleus. The J/?/; 
production is calculated based on 2— >1 kinematics, 

where (M^g) = 2AIj^^ is fixed based on the cc in- 
variant mass distribution predicted by the color singlet 
model. The calculation includes shadowing, taken from 
the nDSg nPDF set, as well as nuclear breakup and the 
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Cronin effect. The nuclear breakup is calculated using a 
parametrization of the dipole cross section fitted to mea- 
surements of the proton structure function at HERA (s^l , 
yielding a breakup cross section that is dependent on 
kinematics of the J/V'- The results from this calcula- 
tion are shown for the 0-100% i?dAu at all rapidities in 
Fig. [S] as the dot-dashed curves. This is a parameter free 
calculation, with no overall normalization or fits to the 
data presented here. The pt shape is in good agreement 
with the data at mid and forward rapidity, but the theory 
shows a greater overall level of suppression than is seen 
in the data. At backward rapidity there is a disagree- 
ment with the shape of the pr distribution. While the 
theory predicts a similar pr shape at all rapidities, the 
data show a much faster rise in RdAu with increasing pT 
at backward rapidity. It is also worth noting, as shown 
in (35| . that this model cannot simultaneously describe 
the rapidity dependence of the PHENIX -R^au and i?cp , 
which is the ratio of RdKn in central collisions to the -R^au 
in peripheral collisions, for J/^ production and therefore 
may not have an accurate description of the geometric 
dependence of the modification. 

A second set of model calculations, performed by Lans- 
berg et al. [13, [s^ , are shown in Figure [9l This model 
uses a Monte-Carlo approach within a Glauber model of 
rf+Au collisions. The J/tj) production is calculated using 
the color singlet model that utilizes 2 — >■ 2 kinematics, 
namely g + g ^ J/ip + g, where the majority of the J ftp 
Pt is balanced by the emission of a hard gluon in the fi- 
nal state, rather than 2 — )■ 1 processes, where the J/4> Pt 
comes entirely from the transverse momentum carried by 
the colliding gluons. The J/^ production is modified in 
d+Au collisions by shadowing effects parametrized using 
various nPDF sets. The calculations shown in Fig. [9] uti- 
lize the nDSg nPDF set. Similar calculations using the 
EKS98 and EPS08 ^ nPDF sets can be found in 
Nuclear breakup of the J/ i/j is taken into account through 
the use of an effective, pT-independent, absorption cross 
section of 4.2 mb. Results using aabs =0, 2.6, and 6 mb 
can also be found in [23|- We have chosen to highlight 
only aabs = 4.2 mb here as it reproduces the rapidity de- 
pendence of the 60-88% RdAu reasonably well [27[ where 
shadowing corrections are expected to be small. The re- 
sults of this calculation, shown in Fig.lHlfor 0-100% -R^au 
at all rapidities, shows reasonable agreement with the 
overall level of modification seen at low pT in the data 
at mid and forward rapidities while the calculation pre- 
dicts a flatter distribution with increasing pT than is seen 
in the data. The shape of the distribution at backward 
rapidity is markedly different than the data. While the 
data rapidly increase to i?dAu « 1 at low pt, the cal- 
culation shows a RdAu that is essentially constant with 
increasing pT ■ 

When comparing the two sets of model calculations in 
Fig. [9] the calculations from Kopeliovich et al. [s^ Issj 
have a different and more pronounced shape when com- 
pared to the calculations from Lansberg, et al [13, . 
Both sets of calculations utilize the nDSg nPDF set, sug- 



gesting a common contribution from shadowing. How- 
ever, the Jf-ip production kinematics are calculated dif- 
ferently, which will lead to some difference in the shad- 
owing contribution. The calculations from Kopeliovich et 
al. include the Cronin effect, which provides a decrease 
in Jf-ip production at low pT and an increase at higher 
Pt, creating an RdAu that exhibits less suppression at 
high Pt than at low pt- The calculations from Lansberg 
et al. do not include the Cronin effect, and therefore the 
Pt shape of iidAu should be dominated by the effect of 
shadowing, and therefore the choice of nPDF set. 

The spatial dependence of the shadowing has been 
taken into account in [27|, where it is assumed that the 
shadowing is proportional to the local density. This as- 
sumption allows for calculation of the RdAu vs pt in dif- 
ferent central! ty bins. The results of the calculation in the 
four PHENIX centrality bins are shown in Figures [TTlfT^ 
and 1131 for backward, mid, and forward rapidity, respec- 
tively. Here we have chosen to include calculations using 
the EKS98 nPDF set along with those using the nDSg 
nPDF set as this will provide a direct comparison be- 
tween the effects due to different nPDF sets, since the 
J/ijj production kinematics and aabs values are identi- 
cal between the two calculations. At mid and forward 
rapidity the calculations are similar to each other and 
show reasonable agreement with the RdAu distributions 
within the current statistical and systematic uncertain- 
ties, although the calculation appears to predict a slightly 
larger average suppression for peripheral collisions at for- 
ward rapidity than is seen in the data. This could be due 
to the value of aabs used at forward rapidity, as the value 
of 4.2 mb was chosen by eye rather than fitted to the 
data, and it may not be independent of y. 

At backward rapidity the calculations are in disagree- 
ment with the data for all but the most peripheral colli- 
sions, where both the calculations and the data show an 
RdAu consistent with 1.0 at all pT- While the calcula- 
tions at backward rapidity using the nDSg nPDF set are 
roughly constant with pT for each centrality, the calcula- 
tions using the EKS98 nPDF set show an enhancement 
in the suppression of RdAu with increasing pT for central 
and midcentral collisions, whereas the data shows the op- 
posite trend. At backward rapidity and low pt (Bjorken 
X « 0.1 for the parton in the Au nucleus) production 
occurs in the anti-shadowing region, while at high pT 
(x ~ 0.3) production begins to move towards what is 
termed the EMC [11] region. In 1986, a suppression of 
the quark distributions within nuclei was discovered in 
the range 0.35 < a; < 0.7 by the European Muon Collab- 
oration (EMC) [3^ in deep inelastic scattering. While 
there is still debate about the source of this suppres- 
sion in the quark distributions, no direct evidence of an 
EMC effect has yet been reported in the gluon distribu- 
tions. Few constraints exist in this region, and there is 
large disagreement in the modification of the gluon den- 
sity between nPDF's. The nDSg nPDF set includes no 
suppression in the EMC region, and only a small anti- 
shadowing effect, while the EKS98 nPDF exhibits a sup- 
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pression in the EMC region similar to that observed in 
the quark distributions, and a larger anti-shadowing ef- 
fect (see [s^l for a comparison of nPDF sets) . The larger 
anti-shadowing combined with the inclusion of an EMC 
effect in the EKS98 nPDF set cause a decrease in the 
calculated RdAu as pr (and correspondingly, x) increases. 
The lack of a strong anti-shadowing effect combined with 
the absence of an EMC effect in the uDSg nPDF causes 
the calculation of RdAu to remain roughly constant with 
increasing px- 

In 39] the authors infer from measurements of T pro- 
duction at RHIC that a strong EMC effect must be 
present to explain the observed modification. Depending 
on the mapping of the J/ip y and pr to x, which is model 
dependent, the high pT data at backward rapidity may 
allow us to probe this region. The large uncertainties 
present in the high px RdAm along with complications 
from competing physics effects in this region, however, 
prevent any strong conclusions from being drawn at this 
time. 
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FIG. 14: (Color Online) J/tp RdAu, as a function of pt mo- 
mentum for midrapidity 0-100% centrality integrated d-|-Au 
collisions. The curves are theoretical calculations from (40l | 
described in the text. 

A third set of model calculations by Sharma and 
Vitev [i^l is compared with the midrapidity 0-100% cen- 
trality integrated RdAu in Fig. [T31 This model describes 
J/ip production using nonrelativistic quantum chromo- 
dynamics (NRQCD). The effect of nuclear shadowing is 
calculated using EKS98 in the EMC region {x > 0.25), 
while for lower values, power suppressed coherent final- 
state scattering leads to a modification of parton x [4l| . 
Initial state energy loss, which accounts for the radiative 
energy loss of the incoming particles through multiple in- 
teractions with the target nucleus is included. This effect 
reduces the energy of the incoming parton, so, to achieve 
the same final-state kinematics the parton must have a 
greater momentum, and therefore a larger value of x. 
This effectively shifts the portion of the gluon distribu- 
tion sampled to higher x. Also included is a calculation 



of the Cronin effect. The solid curve in Fig. [Ml shows the 
full calculation including the Cronin effect. The dashed 
curve in Figure [TJ] is the same calculation without the 
Cronin effect. This comparison gives a direct indica- 
tion of the contribution from the Cronin effect, which 
is evidently over predicted when compared to the data. 
The results presented here will hopefully provide a much 
needed constraint on the Cronin effect at RHIC ener- 
gies. The calculation including the Cronin effect indicate 
an RdAu that decreases at higher px- This is consistent 
with data, however the current statistical and system- 
atic uncertainties make determining the precise trend of 
RdAu difficult at high pT- Better data with a larger px 
coverage is needed to determine the J/ modification at 
higher pT- 



VII. SUMMARY & CONCLUSIONS 

We have measured the J/tp invariant yield and RdAu, 
as a function of px over three rapidity ranges in d-|-Au 
collisions at = 200 GeV using the PHENIX detec- 

tor. These measurements provide a large improvement 
in statistical precision and px reach over the previously 
pubUshed PHENIX d-|-Au resuhs [H, [T^l, and are the 
first measurements of the centrality dependence of the 
J/ px distribution in c?+ Au collisions by PHENIX. The 
A{px) values determined from the data show a marked 
increase with iVcoii that is similar at all rapidities. 

The RdAu vs Px displays similar behavior at mid and 
forward rapidity, showing suppression at low px with a 
gradual increase to a value consistent with 1.0. The RdAu 
at backward rapidity has a different distribution with px, 
showing a more rapid increase from suppression to a value 
of 1.0, and transitioning to RdAu > 1-0 above 2 GeV/c . 
These trends are greater for central collisions, while the 
peripheral collisions show RdAu consistent with 1.0 across 
all rapidities. 

We find an average RdAu for pt > 4 GeV/c of 1.27 ± 
0.06 ± 0.11 at backward rapidity, and an RdAu consis- 
tent with 1.0 at mid and forward rapidity. This implies a 
CNM contribution in A + A collisions that is likely con- 
sistent with 1.0 at high px across all rapidity. This could 
potentially explain the reported increase in Raa with in- 
creasing Px [26[ . However more data and further work 
to understand the propagation of RdAu to Raa is needed 
to confirm this. 

A comparison of the measured -R^Au with three types 
of theoretical calculations was shown. The parameter 
independent dipole model of J/^p production m p + A 
collisions agrees well with the shape of the data at mid 
and forward rapidities, while the shape of the predicted 
Px dependence is different from the data at backward ra- 
pidity. However the suppression is over-predicted at all 
rapidities. The second model uses 2 — >■ 2 J/ip produc- 
tion kinematics coupled with shadowing taken from both 
EKS98 and nDSg nPDF sets as well as an effective ab- 
sorption cross section of 4.2 mb. The calculations with 
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both EKS98 and nDSg show good agreement with the 
data at midrapidity in each centrality bin, as well as the 
centrality integrated case. At forward rapidity the shape 
of the distribution is in reasonable agreement with the 
data, while the overall level of suppression seems to be 
greater in the model calculations than the data. At back- 
ward rapidity, the model calculations using both EKS98 
and nDSG nPDF sets are in strong disagreement with 
the data. At backward rapidity calculations using the 
nDSg nPDF set show a suppression that is constant with 
Pt, while those using the EKS98 nPDF set predict an in- 
crease of suppression with increasing pT- The data show 
the opposite trend. The third model, an NRQCD calcu- 
lation of high PT J I'll} production show a Cronin effect, 
which although generally consistent with the data, is sig- 
nificantly larger than observed in the data, and a sup- 
pression at high pt that cannot be confirmed due to the 
large uncertainties at high pT and the limited pt reach 
of the current data. 

In summary, the data presented here cover a large 
range in x and , providing a further constraint on the 
modification of the gluon distribution in nuclei, as well 
as providing constraints on the size of the Cronin effect 
on J/-0 production at RHIC. 



ACKNOWLEDGMENTS 

We thank the staff of the Collider-Accelerator and 
Physics Departments at Brookhaven National Labora- 
tory and the staff of the other PHENIX participating 
institutions for their vital contributions. We thank Jean- 
Philippe Lansberg, Nicolas Matagne, Boris Kopeliovich, 
Ivan Vitev, and Rishi Sharma for useful discussions and 
theoretical calculations. We acknowledge support from 
the Office of Nuclear Physics in the Office of Science of 
the Department of Energy, the National Science Foun- 
dation, Abilene Christian University Research Council, 
Research Foundation of SUNY, and Dean of the Col- 
lege of Arts and Sciences, Vanderbilt University (U.S. A), 
Ministry of Education, Culture, Sports, Science, and 
Technology and the Japan Society for the Promotion 
of Science (Japan), Conselho National de Desenvolvi- 
mento Cientifico e Tecnologico and Fundagao de Am- 
paro a Pesquisa do Estado de Sao Paulo (Brazil), Nat- 
ural Science Foundation of China (P. R. China), Min- 
istry of Education, Youth and Sports (Czech Repub- 
lic), Centre National de la Recherche Scientifique, Com- 
missariat a I'Energie Atomique, and Institut National 
de Physique Nucleaire et de Physique des Particules 
(France), Ministry of Industry, Science and Tekhnolo- 
gies, Bundesministerium fiir Bildung und Forschung, 
Deutscher Akademischer Austausch Dienst, and Alexan- 
der von Humboldt Stiftung (Germany), Hungarian Na- 
tional Science Fund, OTKA (Hungary), Department of 
Atomic Energy and Department of Science and Technol- 
ogy (India), Israel Science Foundation (Israel), National 
Research Foundation and WCU program of the Ministry 



Education Science and Technology (Korea), Ministry of 
Education and Science, Russian Academy of Sciences, 
Federal Agency of Atomic Energy (Russia), VR and the 
Wallenberg Foundation (Sweden), the U.S. Civilian Re- 
search and Development Foundation for the Independent 
States of the Former Soviet Union, the US-Hungarian 
Fulbright Foundation for Educational Exchange, and the 
US- Israel Binational Science Foundation. 



Appendix A: Details on the calculation of (py) 
1. Fitting The J/ip Invariant Yields 

The J /ip invariant yields as a function of pt were fitted 
with a modified Kaplan function of the form 



fijpr) = po ( 1 - ( ^ 



(Al) 



The data points were compared to the integral of the 
function over the pT bin when calculating the x^. The 
fit results, along with the ratio of the data to the fit are 
shown in Fig. [15] for p+p collisions and Figure [16] for 
0-100% centrality integrated d-l-Au collisions. The fit 
results for each centrality bin are shown in Figs. [TT] [T51 
and ll9l for backward, mid, and forward rapidities, respec- 
tively. 



2. Calculating the Correction Factor k 

To account for the fact that the experimental upper pT 
limits on the J/ip invariant yield distributions vary with 
rapidity and centrality, a correction factor was calculated 
using the fits described in Sec. lA II The ratio 



k = 



(A2) 



was calculated from the fit and applied to the numeri- 
cally calculated {p't)\pt<Pt"'^- The correction factors are 
shown in Table [Vl] and are in all cases small {k < 1.03). 
The uncertainty on k is derived from the fit uncertainty 
by varying the data points within their statistical uncer- 
tainties, refitting, and thereby finding the variation in 
k. 



3. Propagating the Type B uncertainties to {p^) 

When propagating the Type B systematic uncertain- 
ties on the J/i/j invariant yields to the calculated (p^) 
values, the type B uncertainties are assumed to be nor- 
mally distributed. With this assumption we indepen- 
dently sample the uncertainty distribution of the first 
and the last data point of the pt distribution. We then 
assume the Type B uncertainties are linearly correlated 
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FIG. 15: (Color Online) Results of modified Kaplan fits to midrapidity p+p (Left:) and forward rapidity p+p (Right). 




-2.2<y<-1.2, MBdAuf 
X=/NDF=0.64 



1 2 3 4 5 



7 8 9 

P, (GeV/c) 




|y|<0.35, MB dAu 
X2/NDF=1.48 



is ^ 






[ 


1 [ 


1 [ 


1 



s 



2 4 6 



10 12 14 16 

P, (GeV/c) 




1.2<y<2.2, MB dAu' 
X2/NDF=0.83 



1 2 3 4 5 



7 8 9 

P, (GeV/c) 



FIG. 16: (Color Online) Results of modified Kaplan fits to backward rapidity 0-100% d+Au (Left), midrapidity 0-100% d+Au 
(Center), and forward rapidity 0-100% d+Au (Right). 



between these two values. The resulting distribution of 
the values that arises from this procedure gives an 
estimate of the effect of the Type B uncertainties on the 
value of (p^) . 

The Type C systematic uncertainties on the J/ip in- 
variant yields do not affect the calculation of (p^). The 
Type C uncertainties are a global uncertainty, which can- 
cels in the calculation. 
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FIG. 17: (Color Online) Results of modified Kaplan fits to backward rapidity d+Au collisions for each centrality. 
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FIG. 18: (Color Online) Results of modified Kaplan fits to midrapidity d+Au collisions for each centrality. 
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FIG. 19: (Color Online) Results of modified Kaplan fits to forward rapidity d+Au collisions for each centrality. 



System 



TABLE VI: {p^) correction factors, k, for p+p and d+Au collisions. 



y range 



Centrality 



k 



P+P 
P+P 



1.2 < |y| < 2.2 
\y\ < 0.35 



1.006±0.001 
1.018±0.007 



d+Au 
d+Au 
d+Au 



-2.2 <y < -1.2 

\y\ < 0.35 
1.2 < y < 2.2 



0-100% 
0-100% 
0-100% 



1.005±0.001 
1.002±0.001 
1.010±0.002 



d+Au 
d+Au 
d+Au 
d+Au 



-2.2 <y < 
-2.2 <y < 
-2.2 <y < 
-2.2 <y < 



-1.2 
-1.2 
-1.2 
-1.2 



0-20% 
20-40% 
40-60% 
60-88% 



1.005+0.002 
1.002+0.001 
1.006+0.002 
1.008+0.003 



d+Au 
d+Au 
d+Au 
d+Au 



\y\ < 0.35 
\y\ < 0.35 
\y\ < 0.35 
\y\ < 0.35 



0-20% 
20-40% 
40-60% 
60- 



1.001+0.001 
1.002+0.002 
1.024+0.017 
1.020+0.024 



d+Au 
d+Au 
d+Au 
d+Au 



1.2 < y < 2.2 
1.2 <y < 2.2 
1.2 < y < 2.2 
1.2 < y < 2.2 



0-20% 
20-40% 
40-60% 
60-88% 



1.009+0.003 
1.012+0.003 
1.023+0.006 
1.010+0.004 
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Appendix B: Data Tables 



TABLE VII: Data tables for 0-100% centrality-integrated 
RdAu at -2.2 <y< -1.2. 



PT [GcV/c ] 




Type A 


Type B 


TypeC 


0.00 - 


0.25 


0.759 


±0.1 


±0.053 


±0.063 


0.25 - 


0.50 


0.772 


±0.094 


±0.054 


±0.064 


0.50 - 


0.75 


0.853 


±0.092 


±0.059 


±0.071 


0.75 - 


1.00 


0.899 


±0.06 


±0.062 


±0.074 


1.00 - 


1.25 


0.955 


±0.048 


±0.065 


±0.079 


1.25 - 


1.50 


0.934 


±0.048 


±0.064 


±0.077 


1.50 - 


1.75 


0.988 


±0.048 


±0.068 


±0.082 


1.75 - 


2.00 


1.000 


±0.049 


±0.069 


±0.083 


2.00 - 


2.25 


1.043 


±0.054 


±0.072 


±0.086 


2.25 - 


2.50 


1.182 


±0.061 


±0.081 


±0.098 


2.50 - 


2.75 


1.159 


±0.059 


±0.08 


±0.096 


2.75 - 


3.00 


1.161 


±0.067 


±0.08 


±0.096 


3.00 - 


3.25 


1.150 


±0.075 


±0.079 


±0.095 


3.25 - 


3.50 


1.059 


±0.076 


±0.073 


±0.088 


3.50 - 


3.75 


1.234 


±0.101 


±0.085 


±0.102 


3.75 - 


4.00 


1.043 


±0.098 


±0.072 


±0.086 


4.00 - 


4.25 


1.285 


±0.15 


±0.089 


±0.106 


4.25 - 


4.50 


1.133 


±0.152 


±0.078 


±0.094 


4.50 - 


4.75 


1.556 


±0.252 


±0.108 


±0.129 


4.75 - 


5.00 


1.265 


±0.256 


±0.089 


±0.105 


5.00 - 


5.50 


1.186 


±0.303 


±0.083 


±0.098 


5.50 - 


6.00 


1.227 


±0.518 


±0.085 


±0.101 


6.00 - 


7.00 


1.228 


±0.537 


±0.085 


±0.102 


7.{K) - 


8.00 


l.Gi:! 


=1. ():!() 


=().11() 


=0.136 



TABLE IX: Data tables for 0-100% centrality-integrated 
RdAu at 1.2 <y< 2.2. 



PT [GcV/c ] 


RdAu 


Type A 


Type B 


Type C 


0.00 - 


0.25 


0.693 


±0.064 


±0.052 


±0.057 


0.25 - 


0.50 


0.690 


±0.059 


±0.052 


±0.057 


0.50 - 


0.75 


0.664 


±0.048 


±0.05 


±0.055 


0.75 - 


1.00 


0.659 


±0.037 


±0.049 


±0.054 


1.00 - 


1.25 


0.652 


±0.032 


±0.048 


±0.054 


1.25 - 


1.50 


0.671 


±0.032 


±0.05 


±0.055 


1.50 - 


1.75 


0.739 


±0.029 


±0.055 


±0.061 


1.75 - 


2.00 


0.703 


±0.029 


±0.052 


±0.058 


2.00 - 


2.25 


0.732 


±0.031 


±0.054 


±0.061 


2.25 - 


2.50 


0.772 


±0.035 


±0.057 


±0.064 


2.50 - 


2.75 


0.764 


±0.037 


±0.057 


±0.063 


2.75 - 


3.00 


0.821 


±0.043 


±0.061 


±0.068 


3.00 - 


3.25 


0.844 


±0.049 


±0.063 


±0.07 


3.25 - 


3.50 


0.716 


±0.048 


±0.053 


±0.059 


3.50 - 


3.75 


0.765 


±0.056 


±0.057 


±0.063 


3.75 - 


4.00 


0.786 


±0.07 


±0.059 


±0.065 


4.00 - 


4.25 


1.032 


±0.107 


±0.077 


±0.085 


4.25 - 


4.50 


i.UoU 


1 rt 1 on 


_i_n n'7'7 


±U.UoO 


4.50 - 


4.75 


1.118 


±0.166 


±0.084 


±0.092 


4.75 - 


5.00 


1.047 


±0.186 


±0.079 


±0.087 


5.00 - 


5.50 


0.836 


±0.195 


±0.063 


±0.069 


5.50 - 


6.00 


0.984 


±0.298 


±0.074 


±0.081 


6.00 - 


7.00 


1.122 


±0.387 


±0.084 


±0.093 


7.00 - 


8.00 


1.275 


±0.83 


±0.097 


±0.105 



TABLE VIII: Data tables for 0-100% centrality-integrated 
-RdAu at \y\ < 0.35. 



T [GeV/c ] 


RdAu 


Type A 


TypeB 


TypeC 


0.0 - 0.5 


0.78 


±0.065 


±0.11 


±0.061 


0.5 - 1.0 


0.69 


±0.036 


±0.094 


±0.055 


1.0 - 1.5 


0.68 


±0.035 


±0.092 


±0.053 


1.5 - 2.0 


0.78 


±0.048 


±0.1 


±0.061 


2.0 - 2.5 


0.90 


±0.076 


±0.12 


±0.071 


2.5 - 3.0 


0.84 


±0.092 


±0.11 


±0.066 


3.0 - 3.5 


1.00 


±0.14 


±0.13 


±0.078 


3.5 - 4.0 


1.10 


±0.21 


±0.15 


±0.087 


4.0 - 5.0 


1.10 


+0.22 
-0.21 


±0.14 


±0.084 


5.0 - 7.0 


0.91 


±0.21 


+0.11 
-0.13 


±0.072 


7.0 - 9.0 


0.77 


+0.4 
-0.34 


+0.09 
-0.12 


±0.06 
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TABLE X: Data tables for RdAu as a function of pr for 0-20% 
centrality. 







y 




PT [GeV/c ] 




±ype J\ 


j.ype r> 


±ype 


_2 




y < 


— 1 2 


0.00 - 





25 





702 


ibO 121 


±0 049 


±0 063 


2 




y < 


^2 


0.25 - 





50 





760 


itO 104 


±0 053 


±0 068 


2 




2/ < 


2^ 2 


0.50 - 





75 





840 


ibO 105 


±0 058 


±0 075 


_2 


9 <^ 


3/ ^ 


— 1.2 


0.75 - 


1 


00 





894 


±0.076 


±0.062 


±0.08 


2 


9 ^ 


!/ < 


^ 2 


1.00 - 


1 


25 





954 






-I-n n&K 


2 


9 ^ 


y < 


-j^ 2 


1.25 - 


1 


50 


1 


008 




-l-fl riRO 


-I-n no 


2 




y < 


2 


1.50 - 


1 


75 


1 


000 




+n nfiQ 
iu.uoy 




2 




y < 


^2 


1.75 - 


2 


00 


1 


058 


iO 066 


±0 073 


±0 095 


_2 




3/ < 


^ 2 


2.00 - 


2 


25 


1 


058 


iO 062 


±0 073 


±0 095 


_2 


9 <f' 


3/ < 


— 1 2 


2.25 - 


2 


50 


1 


252 


iO 072 


±0 086 


±0 112 


2 


9 ^ 


y < 


-j^ 2 


2.50 - 


2 


75 


1 


240 


-I-n n7R 




-I-n 111 

iU. Ill 


2 


9 ^ 


y < 


"1^ 2 


2.75 - 


3 


00 


1 


249 


-i-c\ nsT 


±0 086 


±0 112 


2 


9 ^ 


y < 


^2 


3.00 - 


3 


25 


1 


230 




-I-n ns 


-I-n 1 1 


2 


9 <-' 


y < 


-^2 


3.25 - 


3 


50 


1 


088 


iO 093 


±0 075 


±0 097 


_2 




y < 


— 1 2 


3.50 - 


3 


75 


1 


406 


iO 13 


±0 097 


±0 126 


-2 


2 < 


y < 


-1.2 


3.75 - 


4 


00 


1 


131 


±0.124 


±0.078 


±0.101 


-2 


2 < 


y < 


-1.2 


4.00 - 


4 


25 


1 


244 


±0.165 


±0.086 


±0.111 


2 


2 < 


y < 


-1.2 


4.25 - 


4 


50 


1 


209 




+n ns/1 


+n 1 ns 


2 


2 < 


y < 


-1.2 


4.50 - 


4 


75 


1 


824 




+n 1 97 


+n 1 fi*? 
loo 


_2 


2 < 


y < 


-1.2 


4.75 - 


5 


00 


1 


508 


iO 329 


±0 106 


±0 135 


_2 


2 < 


y < 


-1.2 


5.00 - 


5 


50 


1 


182 


diO 33 


±0 082 


±0 106 


_2 


2 < 


y < 


-1.2 


5.50 - 


6 


00 


1 


171 


itO 522 


±0 082 


±0 105 


2 


2 < 


y < 


-1.2 


6.00 - 


7 


00 


1 


848 




-I-n 1 on 


-I-n 1 UK 


n 

— z 


2 < 


y < 


-1.2 


7.00 - 


8 


00 


2 


079 


It i.004 


-l-H 1 Ad 


-I-n 1 Qfi 




\y\ < 0.35 


0.0 - 





5 


0.69 


±0.071 


±0.094 


±0.059 




< 0.35 


0.5 - 


1 





0.60 


±0.04 


±0.081 


±0.051 




\y\ < 0.35 


1.0 - 


l 


5 


0.65 


±0.042 


±0.088 


±0.056 




\y\ < 0.35 


1.5 - 


2 





0.73 


±0.054 


±0.097 


±0.062 




|y| < 0.35 


2.0 - 


2 


5 


0.88 


±0.088 


±0.12 


±0.075 




\v\ < 0.35 


2.5 - 


2 





0.81 


±0.11 


±0.11 


±0.07 




Is/I < 0.35 


2.0 - 


3 


5 


1.10 


±0.17 


±0.15 


±0.096 




\y\ < 0.35 


3.5 - 


4 





1.10 


+0.25 
— 0.24 


±0.15 


±0.098 




|y| < 0.35 


4.0 - 


5 





1.10 


+0.26 


±0.15 


±0.095 




\y\ < 0.35 


5.0 - 


7 





0.59 


+0.21 
— 0.19 


+0.073 
— 0.087 


±0.051 




\y\ < 0.35 


7.0 - 


9 





1.10 


+0.61 
-0.54 


+ 0.13 
-0.17 


±0.095 


"1^ 


2 < 


y < 


2.2 


0.00 - 





25 





566 


-I-n r\KA 


-I-n ciAO 


-I-n nf^i 


]^ 


2 < 


y < 


2.2 


0.25 - 





50 





557 


±0 046 


±0 042 


±0 05 


]^ 


2 < 


y < 


2.2 


0.50 - 





75 





557 


±0 034 


±0 042 


±0 05 


]^ 


2 < 


y < 


2.2 


0.75 - 


1 


00 





563 


±0 031 


±0 042 


±0 05 


]^ 


2 < 


y < 


2.2 


1.00 - 


1 


25 





559 


±0 029 


±0 042 


±0 05 




2 < 


y < 


2.2 


1.25 - 


1 


50 





594 


-I-n noQ 


-I-n CiAA 


-I-n ciK'i 




2 < 


y < 


2.2 


1.50 - 


1 


75 





634 


+n n'^^ 


±0 047 


±0 057 


]^ 


2 < 


y < 


2.2 


1.75 - 


2 


00 





587 


±0 029 


±0 044 


±0 052 


]^ 


2 < 


y < 


2.2 


2.00 - 


2 


25 





698 


±0 034 


±0 052 


±0 062 


]^ 


2 < 


y < 


2.2 


2.25 - 


2 


50 





691 


±0.037 


±0.051 


±0.062 




2 < 


y < 


2.2 


2.50 - 


2 


75 





691 


-I-n nyii 


-I-n ciKO 


-I-n nR9 




2 < 


y < 


2.2 


2.75 - 


3 


00 





713 


+n n/i 7 


+n ni^*? 


+n nfi/1 


■|^ 


2 < 


y < 


2.2 


3.00 - 


3 


25 





812 


+n ni^fi 


±0 061 


+n n?*? 


]^ 


2 < 


y < 


2.2 


3.25 - 


3 


50 





631 


±0 052 


_j_Q Q47 


±0 056 


]^ 
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TABLE XI: Data tables for RdAu as a function of pr for 20- 
40% centrality. 
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TABLE XII: Data tables for RdAu as a function of pr for 
40-60% centrality. 
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TABLE XIII: Data tables for RdAu as a function of pr for 
60-88% centrality. 
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